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Multi-particle correlation techniques are frequently used to study jet shapes and yields in hadronic 
and nuclear collisions. To date, a standard assumption applied in such analyses is that the observed 
correlations arise from either jets and associated hard scattering phenomena, or from a background 
component due to combinatorial pairs connected only through whole even correlations. Within 
this assumption of two essentially independent sources, a fundamental problem centers around 
determining the relative contributions of each component. We discuss the methods commonly used 
to establish the background yield in jet correlation analyses, with a full explanation of the absolute 
background normalization technique which establishes the background yield without assumptions 
about the shape of jet correlations. This is especially important in relativistic heavy ion collisions 
where the jet shapes are significantly distorted from the well separated back-to-back di-jets observed 
in proton-proton collisions. 

PACS numbers: 25.75.Bh, 25.75.Gz 



I. INTRODUCTION 

One of the main goals of high energy nuclear physics 
experiments is to produce and study the properties of 
hot and dense nuclear matter in ultra-relativistic heavy 
ion collisions. Probing this matter with strongly inter- 
acting partons produced in hard scattering processes in 
the initial stages of the collision has been a useful tool for 
deducing the properties of the created matter [1, 2, 3, 4]. 
In high energy physics, hard processes are often studied 
using an algorithm to identify and reconstruct jets com- 
posed of partonic fragments to infer properties of the par- 
ent parton. In heavy ion collisions, such measurements 
are extremely difficult due to the large number of total 
particles compared to the multiplicity within a jet. In- 
stead, one-, two- and three-particle observables have been 
used to study hard scattering phenomena. The study of 
particle yields and their correlations, when compared to 
expectations from p+p collisions, provides quantitative 
information about the modification of jet production in 
the presence of hot nuclear matter. 

Two particle correlations have been a particularly im- 
portant tool for describing modification of jets in the 
nuclear environment. Dramatic modification has been 
observed in jet yields and shapes in central heavy ion 
collisions compared to the smaller baseline colliding sys- 
tems. In p+p collisions, two particle azimuthal corre- 
lations have a small angle near-side peak from particles 
which arise from the same jet and an away-side peak at 
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A(f> = 7T from particles which arise from opposing sides 
of a pair of back-to-back di-jets. In central heavy ion 
collisions at ^/sjvat = 20-200 GeV the picture is quali- 
tatively different. The away-side peak is shifted aside to 
A<f> « 2 rad and is much broader than in p+p collisions [5] 
while the near-side peak also contains an elongated struc- 
ture in the longitudinal direction [6, 7]. However, quan- 
tification of these results requires removal of the combi- 
natorial background of particle pairs which are only cor- 
related through whole-event properties. This combinato- 
rial background is large and must be removed in order to 
extract correlations arising from jet processes. The back- 
ground is largest in central events and at lower pt , where 
it can be as much as two orders of magnitude larger than 
the jet signal. 

As the statistical precision of the correlations data 
from relativistic heavy ion collisions increases, better un- 
derstanding of the accuracy and limitations of the exper- 
imental methods used to measure two particle correla- 
tions is needed. To this end we discuss the merits and 
drawbacks associated with methods commonly used to 
extract the jet signal in multi-particle angular correla- 
tions and specifically discuss a method to determine the 
combinatorial background that requires no assumptions 
about the shape of the jet correlations. 



II. TWO PARTICLE ANGULAR 
CORRELATIONS 

The angular correlation technique has been used ex- 
tensively to deduce jet properties in hadronic and nu- 
clear collisions, and is described in detail in several places 



[6, 8, 9]. A brief description of the method is provided 
here. For simplicity, we will focus the following discussion 
on two-particle azimuthal correlations projected onto the 
transverse plane. 

The angular correlation technique statistically exam- 
ines the relationship between particles classified as trig- 
gers (denoted as type A) and partners (denoted as type 
B). When studying jets, trigger particles are typically 
selected to have larger px values than partners and both 
categories usually have pr > 1 GeV/c, though this is 
necessary. 

Due to the back-to-back production of partons by a 
hard scattering process, the distribution of relative az- 
imuthal angles A<p = cj> A — cj> B is expected to peak at 
Acf> = and 7r. However, the measured shape of this dis- 
tribution may be significantly distorted when measured 
in a real detector due to non-uniform angular coverage 
and dead or inefficient areas. The shape of the detector 
acceptance in Acj> can be determined by pairing triggers 
and partners from different events through a process of 
event mixing. A mixed pair contains correlations due to 
the detector acceptance but the physical correlations are 
eliminated. The ratio of the same event pair distribu- 
tion to that of mixed pairs allows pair acceptance effects 
to cancel, leaving a distribution reflecting only physical 
correlations. This ratio is conventionally defined as the 
correlation function, C(A<f>): 
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where n AB is the number of pairs per event for either 
same or mixed pairs and ( ) indicates averaging over 
many events within some centrality selection. [16]. 

Under the assumption that observed correlations 
arise from two independent and separable sources, the 
d( n sVrne) / dA<p distribution consists of AB pairs that are 
correlated with each other by only whole event corre- 
lations (which we call "background" pairs), and those 
where the particles carry additional spatial correlation. 
The source of these additional correlations is generally 
thought of as an association to a particular hard scatter- 
ing, thus these pairs are labelled "jet" pairs [17]. Single 
particles are assumed to be either from jets or some non- 
jet source. Background pairs contain pairs where 1) one 
particle is from a jet and the other is not, 2) pairs where 
both particles are not from jets, and 3) pairs where both 
particles are from jets and not associated with the same 
hard scattering. 
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The two source model can be written as: 

C (Acp) = J (Acj>) + 6 (1 + 2(v£v B ) cos (2A0)) . 
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FIG. 1: (Color online) An illustration of the role of bo in the 
two source model separation of jet and background pairs. 



The quadrupole anisotropy coefficients v\ A,B ^ are taken 
as inputs from independent measurements of type A and 
B particles. Higher order terms of the anisotropy expan- 
sion are smaller and often neglected. In Eq 3, the approx- 
imation is made that the background pair anisotropy is 
equivalent to the product of the single-particle anisotropy 
coefficients: 
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Within the assumption that the azimuthal modulation of 
the background is independent of the jet signal, the fun- 
damental problem of decomposing the correlation func- 
tion into jet and background components amounts to the 
determination of the background level, &o ; &s shown in 
Fig 1. By equating the final background terms in Eq 2 
and 3 and integrating over Acj), using the definition of 
C(Acf>) in Eq 1, the background level can be expressed in 
terms of per-event pair multiplicities as: 
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In practice &o an d thereby (n AB ) has been calculated 
using two approaches: the first involves scaling &o so that 
some portion of the angular distribution is zero, while 
the second uses a calculation to obtain the quantities 
in Eq 5. We follow the conventional nomenclature and 
refer to the former class of methods as the Zero Yield 
At Minimum (ZYAM) [10] approach, and the latter as 
the Absolute Background Subtraction (ABS) [5, 9, 11] 
technique. These methods are discussed in detail in the 
following sections. 

A quantity that is frequently of interest is the per- 
trigger jet pair yield, which describes the conditional jet 
pair multiplicity as a function of relative azimuthal an- 
gle. The term "conditional" refers to the coincidence of a 
trigger-partner pair within some angular region divided 
by the production rate of triggers in the same centrality 
category. It can be shown that the per-trigger yield of 
total pairs (i.e. including jets and background) is related 
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to the correlation function in a simple way: 
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where (n ) is the mean number of triggers per event. 
Since J(A4>) is the fraction of the correlation function 
from jets, the per-trigger yield of the jet pairs can be 
calculated as in: 
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III. THE ABS METHOD 



FIG. 2: (Color online) Schematic of the basic features of the 
£ calculation. 



The absolute background normalization method is 
based on the assumption that the background is com- 
binatorial in nature and that hard scattering results in 
large correlations between the production rates of jet par- 
ticles. The combinatorial background carries only event- 
wise correlations. The background pair production rate 
is given by the product of the single particle production 
rates: (n AB ) = (n A )(n B ). Thus the true combinatorial 
background level in an ideal case is simply: 
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The values of (n A ) and (n B ) are measurable and, in the 
absence of other correlations, an accurate knowledge of 
these quantities is sufficient to determine the background 
level. However, n A and n B are both dependent on the 
event centrality, and this dependence gives rise to a mul- 
tiplicity correlation. More central events typically have 
both larger n A and n B . Because of the correlation be- 
tween n A and n B when events are grouped into centrality 
bins, the number of measured background pairs is larger 
than that expected from Eq 8, (n AB ) > (n A )(n B ). An 
additional correction is needed to account for this effect 
when calculations are made for data selections that span 
a finite centrality range. 



A. Centrality Bias Correction 

We define a scale factor correction, £, for the produc- 
tion rate product that accounts for the covariance effects 
arising from the centrality bias; £ is defined as: 
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The diagram shown in Fig 2 depicts the basic features 
of the procedure to calculate the centrality correction. 
The single- particle production rates, n A and n B , are a 
function of some global property of the collision related 
to particle production, such as the number of nucleons 
participating in the collisions (N part ) or the number of 



binary nucleon-nucleon collisions (N co u). The variations 
of (n A ) and (n B ) are measured over specific intervals in 
these parameters, which are specified by the width of the 
centrality bins used for event mixing. Typically central- 
ity ranges are approximately 5%. The values of N part or 
N co u are based on the results of a Glauber Monte Carlo 
simulation [12]. For the purpose of narration, we discuss 
the calculation using N part , though N co u is of equivalent 
utility, and calculations using this parametrization are 
also made in parallel. In practice, the variation in £ from 
using the two different parametrizations provides a gauge 
of the systematic uncertainty inherent in the method; the 
two parametcrizations bracket the expected scalings of 
hard and soft production that may both contribute to 
the background pair production. Interpolations between 
measurements of {n A } and (n B ) in narrow centrality bins 
are used to estimate the average production of single par- 
ticles at any particular value of N part . 

In a computational algorithm, £ can be calculated by 
throwing an N part value, according to the distribution 
of events within the centrality selection, w glaub , as taken 
from the Glauber Monte Carlo calculation. The average 
number of type A and B particles is determined from 
the interpolated centrality dependence and their prod- 
uct gives the number of combinatorial pairs in the event. 
Events are created in this manner until £ is numerically 
stable. 

The production of A and B particles at a given N part 
is typically modelled with a Poisson distribution. How- 
ever, the details of this functional form do not affect the 
calculation so long as the displacement from the average 
value is independent between triggers and partners. To 
demonstrate this, £ has been calculated for a delta func- 
tion, a step function spanning ±25% the average, and a 
pair of asymmetric delta functions where the production 
at —25% was twice that at +50% the average. The £ 
at 50-60% centrality using these distributions was found 
to be 1.1012(1), 1.1012(2), and 1.1013(2), respectively, 
where the Poisson form gives 1.1010(6). This agreement 
is within the statistical precision of the computational 
tests. Thus, even though Poisson distributions are of- 
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ten used, they are not in general required so long as the 
deviations from average are independent. 

Using the insight that only the average value is rele- 
vant, £ can be calculated equivalently from the Glauber 
distributions and the yield interpolations by summing 
over all N par t values. The correction for centrality se- 
lection becomes a simple matter of finding the event- 
weighted averages of the three functions (n A , n B and 
n A xn B ) for the centrality bin in question. The expres- 
sion for £ is re- written as: 
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where 



values from 2 
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i indexes sequentially over all N part 
A = n A (N* art ). 
For trends of n A and n B that rise (or fall) in concert, 
the value of £ will be always larger than 1. If either n A or 
n B is independent of centrality, the correction is precisely 
1. If for some reason, one trend rises and the other falls, 
£ will be less than one. In practice, the trends of trigger 
and partner production rates with centrality are in the 
same direction and £ is an upward correction on the pro- 
duction rate product. The magnitude of the correction 
depends on how strongly the trends vary across the cen- 
trality bin compared to the yield of the bin. Since parti- 
cle production rises most quickly in peripheral events, the 
magnitude of £ is largest in this region. For the same rea- 
son wider centrality bins require larger corrections than 
more narrow bins. 

We now provide an example of calculating £ by us- 
ing the charged hadron yields published in Ref. [13]. In 
practice, uncorrected yields should be used to determine 
£ in order to properly take into account the multiplic- 
ity dependence of the reconstruction efficiency. Under 
all but extreme cases, the physical centrality dependence 
dominates the value of £ as detector efficiency usually 
varies only slowly with centrality. Therefore the £ trends 
produced here contain the general features of a typical 
calculation. 

Fig 3 shows the Glauber event distributions for vari- 
ous centrality bin selections in both N part and N co u [12]. 
Invariant yields as a function of N part and N co u for part- 
ners, p B 



2.9 GeV/c, and triggers, = 



5.0 GeV/c, are 

shown in Fig 4. The data, {n A } and (n B ), are fit with the 
following two functional forms, chosen for their smooth- 
ness and well-controlled behavior for large N. The in- 
verse tangent, Eq 11, function is referred to as Fit 1 and 
the exponential function, Eq 12, is denoted as Fit 2. 




FIG. 3: (Color online) N par t (top) and N co u (bottom) distri- 
butions from the Glauber Monte Carlo. 



The calculated values of £ from Eq 10 for these trigger- 
partner selections are shown in Table I. For central col- 
lisions £ is a small correction to the background level, 
however since the background level is large compared to 
the jet signal, the effect of including the centrality cor- 
relations on the extracted jet signal is substantial. The 
centrality correction uncertainty is estimated from the 
spread in calculated values using the N part vs. N co u de- 
scription and from using the two choices of functional 
forms. 



TABLE I: £ values for charged hadrons pairs between 5.0 
GeV/c triggers and 2.9 GeV/c partners. 



Centrality 


Npart 


Npart 


NcoU 


N co u 


(%) 


(Fit 1) 


(Fit 2) 


(Fit 1) 


(Fit 2) 


0-10 


1.0041 


1.0048 


1.0041 


1.0050 


10-20 


1.0097 


1.0107 


1.0082 


1.0089 


20-30 


1.0205 


1.0205 


1.0150 


1.0149 


30-40 


1.0369 


1.0353 


1.0246 


1.0236 


40-50 


1.0606 


1.0582 


1.0405 


1.0392 


50-60 


1.1012 


1.1005 


1.0757 


1.0753 


60-70 


1.1825 


1.1873 


1.1604 


1.1639 


70-80 


1.2918 


1.3065 


1.2966 


1.3091 


80-90 


1.3952 


1.4224 


1.4419 


1.4678 



JAB} 



7arctan(/3iV a ) 



(11) 



= 7 (l-e-^") (12) 

where N is either N part or N co u. Sensitivity to the fit 
functional form is assessed by comparison of the resulting 
£ values from use of the two fits. 



B. Other Correlations 

The factorization of pair quantities into singles prod- 
ucts appears often in pair analyses and the degree to 
which the factorizations hold should be examined in each 
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FIG. 4: (Color online) Invariant yield of charged hadrons as a 
function of N par t and N co ii at pr = 2.9 GeV/c and pr = 5.0 
GeV/c. Data are from Ref. [13] and the errors shown are 
statistical and centrality dependent systematic uncertainties. 
Fits are to Eq 11 and Eq 12. 



case. Particle production rates were discussed above. 
Quadrupolc anisotropics in general vary with centrality 
and are subject to an analogous centrality bias correction 
to Eq. 4, though in practice, it has been estimated to be 
a smaller effect than the uncertainties on the anisotropy 
values and neglected [5]. However, particle production 
rates vary more significantly and the correction described 
above is necessary. Additional correlations other than 



centrality variation between triggers and partners could 
require additional corrections of this kind. In practice, 
however, centrality binning correction is the only signifi- 
cant correlation correction to background pair multiplic- 
ities. 

An example of another correlation of particle produc- 
tion rates is the position of the event along the beam 
pipe within a detector with a finite acceptance on this 
axis. Particle production in a symmetric collision system 
is peaked at mid-rapidity and so more pairs will be re- 
constructed when the event is centered with respect to 
the detector than when the event is off-centered. Since 
the overall variation between the two types of events is 
small compared to the kinds of variation seen above and 
can be made smaller with narrow binning, the issue is 
much less of a problem than variations with centrality 
and typically does not necessitate correction. 

Another source could is triggers and partners in the 
background are significantly produced via the decay of 
the same parent particle. In this case, there is a corre- 
lation in the number of background pairs related to the 
parent particle multiplicity within the event. For the 
background in heavy ion collisions, the source is mostly 
pairs from different production centers and is unlikely to 
be strongly influenced by a single decay process, however 
this should be checked with simulations based on the spe- 
cific analysis cuts and detectors used in a measurement. 



C. Pair Cuts 

Pair cuts are necessary in correlation analyses where 
both trigger and partner particles are measured in the 
same detector subsystem. If the trajectories of the trig- 
ger and partner lie sufficiently close together within the 
detector, the deposited signals may interfere construc- 
tively or destructively in their reconstruction depending 
upon the nature of the reconstruction algorithms used. 
The pair efficiency of these effects is extremely difficult to 
model to the precision required in most analyses. There- 
fore pair efficiency is simplified by cutting all pairs that 
fail proximity cuts in the detector subsystems. With 
sufficient detector segmentation, only a small region is 
masked and the percentage of pairs failing the cut will 
be small. 

The influence on pair multiplicity due to rejection of 
pairs with unacceptably small separation is quantified by 
k, the fraction of randomized pairs that survive the pair 
cuts. Some fraction of random pairs will fail the cut due 
to the finite probability for two particles to pass through 
the same region of the detector. This probability can 
be estimated during event mixing. Since the masked re- 
gions represent the same spatial coverage in all events, 
the value of n has no observed dependence on event mul- 
tiplicity. The multiplicity of the background, after taking 
into account losses due to pair cuts, may be calculated 
from the singles distributions and knowledge of the pair 
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(13) 



cut survival probability via [18]: 

{n AB )=tn{n A ){n B ). 

D. Working Equation 



Thus, fully corrected ABS background levels in realis- 
tic scenarios may be calculated in the manner described 
above as: 



/mAB \ 



(14) 



E. Limits of the Method 



The single particle production rate, n , can be writ- 
ten as n A = j A + b A where j A are particles from jets 
and b A are particles from non-jet sources. A similar de- 
composition can be made for type B particles. Using 
this notation, all pairs in the event can be expanded and 
factorized as: 



(n A n B ) 



((j A + b A ) (j B + b B )) 

(J A J B ) + (j A b B ) + (:fb A ) + {b A b B ) 



A-B 



(J A J 



(j A )(b B ) 
+ (j B )(b A ) + (b A )(b B ) 



(15) 



The combinatorial background as estimated in ABS and 
expanded in terms of j and b becomes: 

(n A )(n B ) = (j A + b A )(j B + b B ) 

= (j A )U B ) + (j A )(b B ) 

+ (j B )(b A ) + (b A )(b B ) (16) 

Note that unlike the last three terms, the first term, 
(j )(j B ), is not part of the background. So the ABS 
subtraction produces an extra term beyond the jet sig- 
nal, {j A j B ) , such that: 

(n A n B ) - ^(n A ) {n B ) = (j a j b ) - &(j A ) (j B ) (17) 

For the background subtraction to work without substan- 
tial over subtraction, this extra term is required to be 
small with respect to the jet signal. 



^(j A )(j B )«(j A j B ) 



(18) 



Since hard scattering produces particles at rates deter- 
mined by the characteristics of the scattering itself, like 
momentum transfer, the jet particle production rates for 
A and B particles will be highly correlated with one an- 
other. The presence of jet triggers increases the likeli- 
hood of production of jet partner particles within the 
same event. 



IV. THE ZYAM METHOD 

The Zero Yield At Minimum (ZYAM) methodology 
sets the normalization of the background contribution 
through an assumption that the jet contribution falls to 
zero yield at some point or region in A(f>. 

In addition to the assumption that the jet and back- 
ground correlations arc from essentially independent 
sources and thus separable, the validity of the ZYAM 
method is conditioned upon the existence of (a) one or 
more points with vanishing yield in the actual jet con- 
tribution, and (b) a sufficiently well-sampled correlation 
function that enables a stable and precise determination 
of the minimum value. 

In heavy ion collisions at sufficiently high transverse 

momenta (p^'^ > 4 GeV/c) or inp+p collisions, the jet 
contribution to the correlation function consists of well 
separated near-side and away-side peaks [14]. In these 
cases due to the relative narrowness of the jet peaks, 
there is a broad region over which the background con- 
tribution dominates and can be determined with little 
bias under the ZYAM method. In the case of modified 
shapes in the jet contribution such as those found at in- 
termediate pt in central heavy ion collisions, the ZYAM 
assumption could be broken by jet broadening or other 
modifications creating yield between the locations of the 
near- and away-side peaks. Without independent verifi- 
cation of the ZYAM assumption, the method can poten- 
tially produce unreliable results due to over-subtraction 
in these cases. 

In the simplest, and least reliable, application of the 
ZYAM procedure, the level of the background contribu- 
tion is adjusted (with the harmonic amplitude remaining 
fixed at its measured value) until one measured bin in 
the jet function is zeroed. Clearly, small bins relative 
to features in the jet contribution are needed to limit 
jet contamination of the ZYAM bin. However, division 
of a fixed sample size into smaller and smaller A</> bins 
increases statistical scattering, and hence the degree to 
which the lowest A</> bin is influenced by downward fluc- 
tuations. A slightly more sophisticated method uses the 
average of three neighboring bins in place of a single bin. 
The moving average of neighboring bins attempts to bal- 
ance the effects of jet contamination and statistical fluc- 
tuations (however, depending on the width of the bins 
and the physics of interest, this broader ZYAM region 
could make the assumption of a zero yield region much 
less valid) . The most stable determination of background 
is to fit the correlation function and raise the background 
contribution to touch the fit at the minimum value. As- 
suming a reliable interpolation can be found (which re- 
quires sufficient statistics or outside assumptions), this 
method affords the best reliability against downward sta- 
tistical fluctuations. 

The statistical uncertainty propagated from the ZYAM 
method can be calculated with a simple Monte Carlo al- 
gorithm. The procedure generates simulated correlation 
functions by sampling, bin by bin in A</>, a new point 
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FIG. 5: (Color online) Crosschecks on ZYAM under subtrac- 
tion due to statistical fluctuations for both back-to-back and 
modified jet shapes 

from a Gaussian whose mean and width are given by the 
measured value and its uncertainty. The ZYAM proce- 
dure is then performed on the generated distribution and 
a ZYAM background normalization is extracted. The 
statistical uncertainty of the ZYAM method is thus es- 
timated by the variation of the background level over 
multiple repetitions of this procedure. Such a calcula- 
tion demonstrates that simply taking the statistical un- 
certainty of the one or three points in the first two imple- 
mentations leads to an underestimation of the statistical 
uncertainty, since it does not account for any positional 
shift in Acj) of the ZYAM point. 

Some failures in the ZYAM method can be tested by 
using a known distribution under various levels of sta- 
tistical sampling. There is a strong downward tendency 
in the ZYAM procedure that must be carefully avoided 
in order to extract reliable per-trigger jet pair yield esti- 
mates. The procedure here is similar to that in extracting 
the statistical scatter only it is the average offset from the 
true value that is being examined. Clearly the true value 
in the measured distribution is not known. However, a 
mock distribution similar to the measured distribution 
can be asserted and then tested under statistical sam- 
plings similar to the data. We show the results of two 
such tests of functional forms. One distribution was cho- 
sen to follow a back-to-back shape. The other was given 
an offset away-side peak. Both results appear in Fig 5. 
The choice of functional form does not significantly al- 
ter the resulting under subtractions for the two cases 
tested. The jet to background ratio was not varied in 



these tests. The single-bin implementation of ZYAM is 
shown to be extremely sensitive to catching sizable statis- 
tical fluctuations at low statistics. The three-bin method 
is more robust against fluctuations, but also fails badly. 
The functional fit ZYAM method works best, but is not 
completely robust against failure without the addition 
of sensible constrains on the fits to the correlation func- 
tions. No constraints, such as reasonable jet widths, were 
required in these tests. Unreasonably narrow jet widths 
(less than one A<fi bin) are responsible for the failures in 
these fits at low statistics. 

To summarize, the ZYAM method may fail in two 
cases. The method may over subtract if there is a sig- 
nificant amount of jet yield at the ZYAM point. The 
method also has problems in some implementations with 
under subtraction when applied at low statistics. Given 
the sensitivity of the ZYAM method to low statistics or 
extremely modified jets, it is good practice to confirm 
the results independently with the ABS method in these 
cases. 



V. THE UNDERLYING EVENT 

The ABS background levels in the most peripheral bins 
have been found to lie below the ZYAM background [5]. 
Qualitatively, this is expected because the ZYAM as- 
sumption puts an upper limit on the background level. 
However, it is possible that measurements are also sensi- 
tive to the underlying event, as seen in p+p collisions [15]. 
The underlying event is thought to be composed of ini- 
tial and final state radiation as well as soft parton in- 
teractions besides the one that created the observed jet. 
These multi-parton interactions are not entirely uncor- 
rected with the jet. Furthermore, as the background in 
a small system is the result of very few soft parton in- 
teractions, the multiplicity resulting from a single soft 
interaction to both trigger and partner may become an 
important effect to model. Thus these effects may intro- 
duce additional correlations in the background beyond 
the centrality correlations which are removed by £. 

In large systems, where the background multiplicity 
is almost entirely driven by impact parameter, these 
variations in the combinatorial background play a much 
smaller role in the average background multiplicities. 
Here the difference between ZYAM and the ABS back- 
ground can bracket the uncertainty on the background 
subtraction. The ABS method will underestimate the 
background by not including any underlying event and 
ZYAM will overestimate the background, possibly remov- 
ing some jet signal. In the small systems at higher mo- 
menta, even this extreme in physics assumptions trans- 
lates into a small uncertainty on the extracted condi- 
tional yields. However, small systems at lower momenta 
fair less well and subtractions may produce significant 
uncertainties in the extracted conditional yields. 
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VI. SUMMARY 

The separation of pairs originating in the underlying 
event from those associated with jet production is of great 
importance for two particle correlations in heavy ion col- 
lisions. We have provided the first complete description 
of how to calculate the effect of centrality correlations, 
which enables the combinatorial background to be sub- 
tracted without assumptions about the shape of the jet 
correlations. This method has the additional advantage 
of having small uncertainties, especially in the case of 
statistics limited probes where the uncertainty on the 
ZYAM background will be large. This method could 
trivially be generalized to three or more particle correla- 
tions. Proper, consistent treatment of the background is 
crucial to the quantitative extraction of properties of the 
hot dense nuclear matter and its geometry from multi- 
particle correlation measurements. 

We have described a method for subtraction of the 
combinatorial background under the assumption that 
there are two independent sources of pairs. Reality, how- 
ever is likely more complicated. The discussion of the 
underlying event in Section V contains one example of 
additional correlations which could be important. Cur- 



rently, the two source assumption has proven useful in the 
interpretation of two particle correlation measurements. 
More precise measurements and more complicated exper- 
imental questions could necessitate a re-evaluation of the 
two source model. In the meantime, it is recommended 
that publications include full correlation functions in or- 
der to make available the data containing no assumptions 
about jet-background decomposition should these consid- 
erations prove necessary. 
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